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Abstract 



Formulas are developed for use in lattice studies of charged hadron polar- 
izabilities. In particular, the valence quark different-flavor component of the 
charged pion polarizability is examined on a 16 3 x 24 lattice at f3 = 6.0 using 
Wilson fermions. Using the elastic limit results of Part I of this series, it is 
concluded that this represents a small negative component of the full charged 
polarizabilty. 
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I. INTRODUCTION 



In the previous paper of this series, referred to here as Part I, the elastic limit of pions and 
rho mesons and were examined and the results were compared to lattice simulations using 
three-point function techniques. The elastic "baseline" set in Part I makes the measurement 
of nonelastic properties of these mesons possible. 

Neutral hadron polarizabilities have been examined previously using the techniques of 
lattice QCD [J]]. The authors there used an external field method to study the effect of a 
uniform electric field on the correlation functions of mesons and baryons. However, as of yet 
there have been no investigations of charged hadron polarizabilities using the lattice tech- 
niques. This is unfortunate since there is now a growing body of experimental information 
on this subject. Of course, external field methods would be very difficult or impossible to 
use in the case of charged hadrons because charged particles accelerate in an electric field. 
The definition of charged hadron polarizability takes the form of the response of the particle 
to an applied oscillating electromagnetic field. For this purpose a single photon will do and 
the effect being considered is then just Compton scattering. 

The subject of charged pion polarizability will be considered from the point of view of 
an effective relativistic field theory. Comparing the general and phenomenological forms of 
the Compton scattering amplitude to second order in the photon momentum will give us the 
form of the polarizability coefficient. Next, the method for extracting polarizability on the 
lattice is explained. It is pointed out that the full pion polarizabilty requires the separate 
evaluation of five types of connected and disconnected fermion diagrams. In the fourth 
section the details of of our numerical simulation of the valence different-flavor component 
of the charged pion polarizabilty will be given and the results from Part I will be extensively 
used. It will be established that the resulting contribution is a small negative component of 
the full charged polarizabilty. 
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II. THEORY 



Charged hadron polarizability has been examined previously in Ref. |1 using the methods 
of charge overlap in the context of lattice QCD simulations. However, the methods of Ref. |2j 
are coordinate-space based; the present techniques || are momentum-space based. First, let 
us review the definition of charged pion polarizability as measured in Compton scattering. 
The definition of charged polarizability in terms of relativistic matrix elements in the context 
of an effective meson theory will then be established. 

A. Review 

The general form of the charged pion Compton scattering amplitude, 

T^ = iJ d'xe^inip^TU^UO))^)), (1) 
where 'T' denotes time ordering and 

Mx) = q*3l(?) + qd jfa), (2) 

(q u — = — |) is an electromagnetic current density for up and down quarks, to second 
order in the photon momentum is given by 

j 2E 2E rp = gUgi + fci,pi)T„(p2,P2 + h) _ T^(p 2 ,p 2 - fci)T„(pi - k 2 ,pi) 
* ^ (pi + ki) 2 — m 2 (pi — k 2 ) 2 — m 2 

+2g tlu + A{k\g liV - k^k lfl + kjg^ u - k 2fl k 2fl ) + B(k 1 - fag^ - k 2fl k lv ) + Ct^. (3) 

In this context A, B and C are simply numerical coefficients. A noncovariant continuum 
state normalization, 

E/(^Kp))Kp)I = i> (4) 

is being used which is responsible for the square root factor involving pion energies appearing 
on the left hand side of Eq.([3]). The pion electromagnetic vertex (g M = p'^ — p^), 



Up',p) = {p', + p,)fm 2 ) + - fAq 2 )), (5) 
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has been written so as to obey 

%T^(p\p)=p 2 -p 2 , (6) 

for off-shell pions, which is needed in order to satisfy the scalar electromagnetic Ward identity 
||. To the fourth order in momentum (needed later), 

Kfi. + <£>S. ( 7 ) 
In the above m is the pion mass and the tensor t^ v is given by 

V = ■ k 2)QvQv + (Q ■ h)(Q ■ h)g^ - (Q ■ k 2 )Q^k lv - (Q ■ k x )Q v k 2ll . (8) 

(Qfi = Pifi + P2fi-) The Compton amplitude is conserved {k^T^ = k^T^ = 0), which can 
be shown to result in A = <r 3 2> to this order. The first two terms on the right of Eq.@ 
represent elastic scattering from the pion while the others are contact terms. (A metric 
with g 00 — 1, g n = — 1, (i = 1, 2, 3) is being used.) In order to define electric and magnetic 
polarizabilities, let us examine the Compton scattering amplitude ae^T^e^*, where a is the 
fine structure constant and e\ and e 2 * are the initial and final photon polarization vectors, 
respectively. In the lab frame where the inital pion is at rest: 

h ■ ei = k 2 ■ e* = 0, (9) 
Pi -ei= Pi •€* = (). (10) 

Choose ei and e 2 to have only spatial components. Defining k9 12 , = ^(1,2), m the nonrela- 
tivistic limit (|A;(i,2)|, ^(1,2) << m ), one can show for off-shell photons that 

aeZT^e™ = h ■ e* 2 [ (1 + — - — (kf + kl)) + a E uJiu 2 } + P M (h x k x ) ■ (e* x k 2 ), (11) 

mo 

where 

a E = -a{— + 2mC), (12) 
2m 

Pm = a—, (13) 

are by definition the electric and magnetic polarizabilities, respectively. An explicit form for 
a E will now be exhibited and an explanation of how this coefficient can be measured in a 
lattice simulation will be given. 



B. Derivation 



First of all, let us use the continuum to lattice corespondences, 



Hp)) — * V 1/2 \n{p) >, (14) 
Ux) — > 3 -M, (15) 



P POO 

/ d 4 x — > Va / dt£ (16) 

where = N s a 3 , N s being the number of spatial sites in the lattice and 'a' being the spacing, 
to rewrite Eq.(P as 

T, v = zN'a J dt^e^ < p 2 |T(^(x)^(0))| Pl >, (17) 

X 

Notice that time remains continuous (and is dimensionless) in this formalism and that a pos- 
sible finite renormalization factor has been neglected in Eq.([L5|) since the exactly conserved 
lattice charge density will be used here. 

Next, a point about the amplitude being calculated needs to be made. On the lattice, 
the appropriate object to consider is actually a normal ordered product of currents: 

Tg' = iN 2 s a J dtJ2e lk2X < 7r(pa)| : #(*)#(0) : |tt( Pi ) > . (18) 

X 

This is because the vacuum amplitude lattice Green function used (see Eq . (|36|) below) 
contains a disconnected part when p2 = pi which must be removed by hand. The form 
of the disconnected amplitude being subtracted in Eq . fll8|) will later be exhibited explicitly. 

The methods of charge overlap will be used to calculate the necessary amplitudes. Thus, 
we consider the four-point function corresponding to Too with the kinematical conditions: 

uji = u 2 = 0; pio = P20 = m; \ki\ = \k 2 \ « m. (19) 

The relation between normal-ordered and time-ordered currents is 

: #(*)#(0) : = T(#(s)#(0)) - < 0|T(^(x)^(0))|0 > . (20) 

Using lattice completeness, 



J2\n(p n )><n(p n )\ = l, (21) 

n,p n 

as well as space-time translations, 

J» = e**#(0) e -*-«, (22) 

the space sum and time integral in Eq. (|T8D now results in 
T e//, _ nAT 2^ \<A0)\p L (0)Hk)>\ 2 N s K0b*i>f ^ 

oo I -u^ - ^ ^ {En _ m) ^ • l^J 

Contact between Eq.(|23"D and Eq.([3D above now needs to be made. Separate off the 
elastic part in Eq. (|23|) , 

<^l^(0)l^)>^^^), (25) 

where F^q 2 ) is the usual continuum pion form factor. Using Eq.(^), one can show after 
some algebra that 

4m r 14 9l -* 9r <r 2 >l 9 1 A , . , 

T n e ' = — + { m<r 2 >} + k 2 { + -m < r 2 > 2 +—m < r 4 >}, 26 

00 k 2 m 3 1 1 3m 9 15 1 V ; 

up to second order in the photon momentum. For the nonelastic part one may show at low 
momentum (\k\ << m) for vector or axial vector states that 

< 7r(0)|/(0)Mfc) >= 6 ° W ^" m " , (27) 

< |/(0)KA) >= ^)^'\ (28) 



which serves to define the constants d nn and do n , and e°(/c) is the 0th component of a 
polarization vector for the state \n(k) > with mass m n . The d nn , d n0 factors are defined in 
these expressions so that they are independent of the lattice size and spacing. Interestingly, 
vectors and axial vectors are the only types of intermediate states which may contribute to 
the polarizability. See the Appendix for further details. If we choose a helicity basis where 
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the e Al (/c, A) with A = +, — are purely spatial and the positive-z axis is along k, then only 
the longitudinal polarization state, A = 0, contributes in Eqs.fl2"7|) and (|2"g|), 

e°(^0) = H (29) 
m n 

and the expansion to second order in photon momentum for these kinematics then yields 

// l ^2 M n\ 2 \don\ 2 

T 00 {"kinematics" = ^00 + ^k ( }^ — — — ■ — V — ), (30) 



where Tqq is given by Eq. fl26|) . 

The phenomenological expression, Eq. (|3"0"D , is to be compared with the \i — v — compo- 
nent of the general Compton amplitude, Eq.(^), also to order k with the same kinematics. 
One may show that 

T (pi + h,Pi)T {p2,P2 + h) T (p 2 ,p 2 - fci)T (pi - k 2 ,Pi) 



2mT ( 



00 1 "kinematics" 



2 



(pi + ki) 2 — m? (pi — k 2 ) 2 — m 

+ 2 - 2<r2> k 2 - Bk 2 - Am 2 k 2 C. (31) 



For these kinematics, one has 

1 1 1 



(32) 



(Pi + h) 2 - m 2 (pi - k 2 ) 2 - m 2 k 2 ' 
and also, for example, 

2 4 

U Pl + h, Pl ) = 2m(l - ^j^q 2 + ^T^ 4 )- (33) 
When the pieces are put together, one obtains 

I 4m ,14 2ir2r <r2> l 22 
Tool "kinematics" = + { ~m < T >} + k { h -771 < T > 

k 2 77i 3 3m 9 

+ — m < r 4 > -(— + 2mC)}. (34) 
15 2m 



As the last step in the derivation, compare Eq. (|30T ) and Eq.(p4|) and identify 

a E = -a(JL + 2mC) = 2a:{£ J^- - ^£ ^}. (35) 

2m n^Tr \ h n -m) n hi., 
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This formula relates to properties of relativistic matrix elements, which was the goal. 
Of course, the d nn and ofon factors are nothing but appropriate derivatives of the matrix 



of this formula worthy of mention. First, the explicit appearance of a volume factor makes 
it clear that neither term above is finite, but the whole expression must be regulated (as 
on the lattice) to extract a finite answer. Second, it does not have the classical charge 
radius term, a< r 2 >/3m. This is because it comes from comparing to an effective theory 
of point mesons. (Even the pion is viewed as a point object in this case; it's "size" simply 
comes from additional local interactions with photons.) Third, it explicitly includes the 
vacuum polarization contribution first explained by Terent'ev in Ref. |J. The negative 
sign in front of this term makes it abundantly clear why in relativistic field theories the 
charged particle's polarizability can be positive or negative. The polarizability of the pion, 
no longer considered the groundstate as in nonrelativistic quantum mechanics, is actually 
defined relative to the vacuum, a result remininiscent of the Casimir effect. Finally, even 
though we are only keeping terms to order k 2 , it is clear that the above is an exact formula 
for the polarizability, which is simply the numerical coefficient of the second order energy 
term in Eq.flTT|). 



Let us concentrate now on the lattice amplitudes necessary to directly extract the charged 
pion polarizability. Some comments about indirect lattice measurements of this quantity will 
be made later. 

For these purposes, one forms the usual charge overlap function, but with normal ordering 
of the charge densities instead of time ordering: 



where £3 > (£2,^1) > 0, z = (0,z) and r = xi — x±. Dimensionless times and distances 




There are several other aspects 



III. LATTICE MEASUREMENT 




(36) 



S 



are being used in this section. The formulas are now in discrete Euclidean space but one 
imagines the time axis is still continuous. The time limits t 3 » ti ;2 » 1 yield 

P w (r;t 3 ,t2,fi) ^< tt(0)| : p L {x 2 )p L {x l ) : |vr(0) >, (37) 

where 

< tt(0)| : p L (x 2 )p L ( Xl ) : |tt(0) > = < tt(0) \T(p L (r)p L 0)) |tt(0) > 

-<0|T(p L (r)p L (0))|0>. (38) 

Consider now the finite Fourier transform of P n (r; t 3 ,t 2 ,ti), 

Q{q; t 3 , t 2 , h) = N S J2 e'^P^f] t 3 , t 2 , h). (39) 

In the same time limit: 

Q(g;t 3 ,*2,ti) — ^E e "^{< 7r(0)\T(p L (r)p L 0))\7r(0) >-< 0\T(p L (r)p L (0))\0 >}. (40) 

r 

Assuming t 2 — t\ > and inserting a complete set of states now gives: 

g(g;t 3 ,*2,*i) — > I <n(0)\p L (0)\n(q) > | V(^-m.a)(t 2 - tl ) 

n 

- E I < °IP L (°) W<?) > pe"^' 2 -' 1 )}. (41) 

n 

Separate off the elastic term (n = 7r), 

Q el (q;t 2 -h) = N*\ < 7r(0)|p L (0)|7r(g) > |2 e -(^«-«»-«)(* a -ti) ) ( 42 ) 



above. Then, using Eqs.(^), (28) and ( p9[) at low momentum transfer results in the con- 
struction, 



/oo 
d(h - h)[Q(q; t 3 , t 2 , h) - Q el (q; t 2 - t,)} 
-oo 



= 2«fc 2 (E - V ^> ^ = ^ (43) 

(A; = <f/a) when the relative time integral between the two charge densities is explicitly 
done. (One gets a doubling of the result when the other time relationship, t 2 — t\ < 0, is 
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assumed.) Of course on a finite lattice the time integrals do not really extend to ±00 because 
of the fixed time locations of the particle sources, but instead to time limits such that the 
elastic limit is effectively established. Eq.(fl^) now gives us the means of measuring the 
Compton scattering pion polarizability coefficient. On the lattice one may take a numerical 
"derivative" of the left hand side of Eq.p3|) with respect to q 2 by simply dividing it by <f 2 at 
the lowest lattice momentum. (The relative error in this procedure is of order q 2 .) 

It is crucial to work in momentum space to extract the polarizability. In position space 
one would look at (image corrections are assumed), 

R 2 (t 3 ,t 2 ,tx) = Y,? 2p ^hMM), (44) 

r 

which does not project onto good momentum when a complete set of states is inserted 0. 
Thus, the limitation to small momenta, necessary to use Eqs. fl27|) and (|28|) above, can not 
be enforced. 

Although Eq.(^3|) is a simple formula, it's direct evaluation on the lattice via four-point 
functions which arise from the field contractions in Eq. ( |5o| ) is not easy. There are five types 
of diagrams, which may be classified by fermion line topology, as shown in Fig.l. In the 
case of Figs. 1(a) and (b), the photons interact directly with the valence quarks. Fig. 1(a) 
will be termed the same-flavor valence part and Fig. 1(b) will be called the different-flavor 
valence part. Figs. 1(c), (d) and (e) represent contributions from disconnected fermion lines 
in which arbitrary numbers of gluon lines connect the various pieces together. These gluon 
connections are enforced in the correlations between the pion propagator and the various 
fermion loops in the lattice Monte Carlo. (When isolated groups of fermion lines arise, 
"connected part" means "statistically correlated part" on the lattice.) Fig. 1(c) represents 
the contribution from correlations between a valence quark and a self-contracted quark 
loop. (Flavors other than u and d can of course contribute to the quark loops.) Fig. 1(d) has 
both currents, either same-flavor or different-flavor, self-contracted. In the case of Fig. 1(e), 
one has a complete quark loop formed. It is the amplitude where the pion propagator 
is completely disconnected from the quark loops in (d) and (e) (when p 2 = Pi) which 

10 



needs explicit subtraction and is responsible for the negative contribution on the right hand 
side of Eq. (|35|) . Fig. 1(c) also has a disconnected part, however, there is no need to do a 
subtraction since it can be shown that the vacuum expectation value of the exactly conserved 
lattice charge density vanishes configuration by configuration for color SU(2), and in the 
configuration average for SU(3). 

The lattice formula for has now been derived and the agenda for evaluating the 
charged pion polarizability directly on the lattice has been set. In the next section we 
will examine one particular type of diagram and show that it actually makes only a small 
contribution to «b. 



IV. LATTICE SIMULATION 

As pointed out immediately above, the direct calculation of the charged pion polariz- 
ability involves many diagrams which must be separately evaluated. In this section the 
different-flavor valence piece, Fig. 1(a), will be considered. 

The full different-flavor piece involves Figs. 1(a), (c) and (d) and can be evaluated on the 
lattice at low lattice momentum, q, by 



u,d <-^ a 
Op = 



/oo 
d(t 2 - h) [Q u > d (q; h, t 2 , h) - Q el ^ d (q; t 2 - h)}, (45) 

where Q u,d {q]t s ,t 2 ,ti) and Q el,u,d (q;t 2 — £1) are identified from the Fourier transform of 

P n (r, t 3 , t 2 , h ) = 2q u q d — < |Q > , (46) 

as above. 

We now specialize to evaluating the valence u, d contribution from Fig. 1(a). We refer the 
reader to Part I for a discussion of the details of the numerical simulation, which was carried 
out on 20 quenched 16 3 x 24 lattices with Wilson fermions at /3 — 6.0. The appropriate 
correlation functions for three quark masses are displayed in Fig. 5 of Part I. According 
to Eq.(H3|), we now need to evalute the area under these functions and subtract the elastic 
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contribution. The results of this procedure are presented in Table I, where the errors specified 
are purely statistical. In arriving at these results, the numerical value aa? = 1.06 x 10 _44 cm 3 
has been used, where a is the fine structure constant. This value is inferred from the inverse 
lattice spacing, a -1 = 1740 MeV, which is taken from comparing the chirally extrapolated 
lattice nucleon mass with the experimental value ||. 

When the results in Table I are extrapolated linearly in quark mass, as defined in Part 
I, to the chiral limit at n cr = 0.1564, we obtain (a u /) valence = -.12(11) x 10" 4 fm 3 . 

There are various systematic uncertainties in our final numerical result. One systematic 
uncertainty is associated with the numerical evaluation of the discrete nonelastic time cor- 
relation integral. This discrete time integral was done using both the trapezoidal rule and 
Simpson's method [0]. The differences for all three k values were less than 2%. There is an 
additional systematic uncertainty associated with the numerical momentum derivative. As 
pointed out above, this is expected to be of order <f 2 , or about 10%, the same as other finite 
lattice spacing errors. This was examined numerically by repeating the simulation (with the 
same fits) on 10 configurations, also at (3 = 6.0, on a larger 20 3 x 30 lattice at k — .152. The 
result on the larger lattices is —.33(14) x 10~ 4 fm 3 , which unfortunately has large statistical 
errors, but is consistent with the value in Table I. In addition to these uncertainties, there 
is also the systematic uncertainty associated with the setting of the lattice scale. Since the 
pion polarizability scales as a 3 , it is rather sensitive to a change in the scale. In order to 
remain consistent with the low energy lattice results for nucleon mass, pion decay constant 
and proton electric form factor in Ref. ||, a roughly 10% change in the present lattice scale 
could be tolerated. This means a systematic scale uncertainty of approximately 33% for «£. 

V. SUMMARY AND CONCLUSIONS 

The establishment of the charge overlap elastic limit in Part I allows the possibility of 
obtaining nonelastic hadron properties directly from lattice QCD. In particular, using the 
lattice formulas developed, we have see that the charged pion polarizability can be measured 
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on the lattice from low momentum charge overlap correlation functions. 

We have obtained (a^ d ) valence = —.12(11) x 10~ 4 fm 3 in the chiral limit from the quark 
diagram in Fig. 1(a). This is to be compared with the experimental result || 

a e * p = 2.2(16) x l(T 4 fm 3 , 

and the result from second order chiral pertubation theory fllCf] , 

a chvral = x 10 -4 fm 3 

Thus, the present results establish that the valence quark different-flavor diagram, Fig. 1(a) 
contributes only a small, negative component to the full charged pion polarizability. We 
must therefore look to the other diagrams in Fig.l for the bulk of this effect. 

It is important to note that in addition to the present, direct measurement of the charged 
pion polarizability, there is also the Das, Mathur, Okubo (DMO) current algebra sum rule 
fTT| for otE- Although this formula is not exact, the measurements required for it's evaluation 
on the lattice are much simpler than the direct method since it involves only mesonic two- 
point functions. A lattice evaluation of the so-called intrinsic part of the charged pion 
polarizability in the DMO sum rule should be quite interesting and is currently underway. 

Finally, it was noted above that the lattice measurement of ag scales like a 3 and is 
therefore quite sensitive to the lattice scale. This opens up the possibility of using charged 
pion polarizability as a place to accurately test QCD, or alternately, to set the lattice scale, 
once the lattice and experimental results are both improved. 
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Appendix 

In this Appendix it is established that the only intermediate states which may contribute 
in Eq.(^) to the charged pion polarizability are vectors and axial vectors. 
First, let us dispose of the spin zero case. We have (n ^ n) 

11 W fnW l \\ PaGljq 2 ) + PnaG 2 {q 2 ) , , 

{n(p)\ja{0)\n(pn)) = (2E 2E n y/ 2 ' ^ ' 

where as usual q 2 = (p — p n ) 2 . Charge conservation requires the right hand side of Eq.(f4"T|) 
dotted into p% — p^ to be zero, so that 

G 2 (q 2 ) = G 1 (q 2 )( m l~ p - Pn ). (48) 
m 2 n -p-p n 

Therefore (considering the pion at zero spatial momentum) 

M0)M0)Kp„)) = (^gv,(- h- (49) 

At low momentum, \p n \ « m, we then find 

(7r(0)|io(0)|n(p n )) oc G 1 ((m - m n f)f n . (50) 

Therefore the absolute square of this matrix element goes at least like p£ at low momentum, 
assuming Gi((m — m n ) 2 ) is nonsingular, and can not contribute to the pion polarizability 
term in Tqq^, which is proportional to p 2 . 

We now consider general nonzero spin, s. Some fundamental properties of polarization 
tensors are: (A s = — s, , , 0, , , s): 

e^(p,\s)e; i ...Jp,K) = Sx s>K , (51) 
W^-'^^-'^AsHO, (52) 
P #te e»-"-"*( P ,A,) = 0. (53) 

Adopting a helicity basis, the Lorentz structure for spin s must be of the form 

(n(p)\ 3a (0)\n(s,\s,Pn)) = ( 54 ) 
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where it is assumed that the tensor A lll ... llsa {p,p n ) is nonsingular. 
First, we will need the following result: 

e -°(p,A s ) oc \p\ s , A s = 0, 
e°-°(p,A a ) = 0, A s ^0. 



(55) 



This can be established from the higher spin considerations in Ref. |12|. This reference 
shows that the polarization vectors for higher spin may be constructed from those of spin 1 
via 



E ^(O^-^a 



Ai=-1 



(56) 



where 



(2s)! 



1/2 



(57) 



{s + \ s )\{s-\ s )\ 

Using a purely timelike and recalling that e°(p, Ai) ^ only when Ai = 0, it is easy to 
see by matching coefficients of £ + and £_ on both sides of Eq.|56| that 



> . 



(58) 



e°-°(p,A,) = 0, A s ^0. 

Now using e°(p, 0) cx \p\ from Eq.(|29|) establishes Eq.(|55|) above. There are many other ways 
of proving this result. 

Working in the rest frame of the pion and considering only jo, the matrix element in 
Eq.(|51D is now seen to be proportional to e°"'°(p, A s ): 



(7r(0)|i (0)|n( S ,A s ,p n )) oce°'-°(p n ,A s ). 



(59) 



This comes about through repeated use of Eqs.(^) and (|53f) . For example, consider the 
s = 2 case. Explicitly, one has 



{7T(6)\j (0)\n(s,X s , Pn )) 



(2m2E n y/ 



00 



^OioO + e ^000 



(60) 
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The remaining spatial 0(3) symmetry of the matrix element requires that A^qq oc {p n )i 1 , 
A 0i20 oc (p n ) i2 , and A ili20 oc {p n )iAPn)n ox A ili20 oc g ili2 . Through the use of Eqs.(0) and 
(|53|) each term in Eq.(|60|) is seen to be proportional to e 00 (p n , A2). This in turn means that 
the matrix element couples only to A2 = and behaves like p% (or possibly higher power) 
at low momentum from Eq.fl55p. This argument can be extended to higher spin. Thus, 
s = l,Ai = (axial) is the only possible contributor to the charged pion polarizability 
through this matrix element. 

Although only the matrix element (7r(0)\jo(0)\n(s, X s ,p n )) has been discussed, similar 
conclusions hold for (0\j Q (0)\n(s, X s ,p n )); namely, that the state n(s,X s ,p n ) contributes to 
charged pion polarizability only when s = 1, Ai = (vector). 
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Table 



TABLE I. Valence contribution to the dfferent-navor charged pion polarizability as a 
function of k. 



k as (fm 3 ) 

0.148 -0.29(6) x 10~ 4 

0.152 -0.23(6) x 10~ 4 

0.154 -0.15(10) x 10- 4 
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Figure Caption 



1. Graphical representation of the four-point functions contributing to charged pion po- 
larization. The zero-momentum pion interpolating fields are represented by horizontal 
bars, the time axis runs horizontally and the vertical direction represents one space di- 
mension, (a) and (b) represent different-flavor and same-flavor valence contributions, 
respectively. Fig. 1(c) is the contribution from a single isolated quark loop, while (d) is 
the two quark loop contribution, (e) arises from a same-flavor current- current contrac- 
tion. It is assumed that gluon lines connect the various fermion lines together in each 
diagram. Diagrams (d) and (e) have a disconnected part which must be subtracted. 
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(a) 




(e) 

Fig.l 



